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Abstract 

We describe a method that allows for a practical application of the theory of pseudo- 
Hermitian operators to PT-symmetric systems defined on a complex contour. We apply 
this method to study the Hamiltonians H = p 2 + x 2 (ix) u with v £ (— 2,oo) that are 
defined along the corresponding anti-Stokes lines. In particular, we reveal the intrinsic 
non-Hermiticity of H for the cases that v is an even integer, so that H = p 2 ±x 2+u , and 
give a proof of the discreteness of the spectrum of H for all v G (—2, oo). Furthermore, 
we study the consequences of defining a square-well Hamiltonian on a wedge-shaped 
complex contour. This yields a PT-symmetric system with a finite number of real 
eigenvalues. We present a comprehensive analysis of this system within the framework 
of pseudo-Hermitian quantum mechanics. We also outline a direct pseudo-Hermitian 
treatment of PT-symmetric systems defined on a complex contour which clarifies the 
underlying mathematical structure of the formulation of PT-symmetric quantum me- 
chanics based on the charge-conjugation operator. Our results provide a conclusive 
evidence that pseudo-Hermitian quantum mechanics provides a complete description 
of general PT-symmetric systems regardless of whether they are defined along the real 
line or a complex contour. 



PACS number: 03.65.-w 



1 Introduction 



The notion of a pseudo-Hermitian operator as outlined in pi |2j Ej provides a general frame- 
work for understanding the intriguing mathematical properties of PT-symmetric Hamilto- 
nians [U Ej. 1 It involves an underlying Hilbert space TC in which the operator acts. For 
PT-symmetric Hamiltonians defined on the real line, TC is the familiar space of square- 
integrable functions. For the PT-symmetric Hamiltonians H defined on a complex contour 
and having a discrete spectrum, TC is the Hilbert space obtained by Cauchy-completing the 
span of the eigenfunctions of H with respect to an arbitrarily chosen positive-definite inner 
product fTJ QUI E] • The implicit nature of this construction makes a direct application of 
the theory of pseudo-Hermitian operators for these Hamiltonians intractable. This forms the 
basis of the view that this theory is incapable of dealing with PT-symmetric Hamiltonians 
defined on a complex contour. The purpose of this article is to show that indeed the opposite 
is true. This is done by an explicit construction that allows for the description of the same 
system using the information given on the real axis. It reveals the implicit non-Hermiticity 
of the apparently Hermitian PT-symmetric Hamiltonians, such as p 2 — x A , that are defined 
along an appropriate complex contour |U|3]. Furthermore, it leads to a previously unnoticed 
connection between the spectral properties of the PT-symmetric Hamiltonians of the form 

H = p 2 + x 2 (ix) v , 

(defined on an appropriate contour) and those of the Hamiltonians of the form 

H = p + \x\ , 

(which are obtained by requiring the eigenfunctions to belong to T 2 (R) and satisfy certain 

boundary conditions at x — 0.) An important advantage of a pseudo-Hermitian description 

of PT-symmetric systems defined on a complex contour is that it offers a prescription for 

computing the physical observables [TJJ EH E] of these theories. 

In the remainder of this section we include a brief review of the relevant aspects of the 

theory of pseudo-Hermitian operators. For clarity of the presentation we will only consider 

lr The term "pseudo-Hermitian" has been in use within the context of indefinite- metric quantum theories 
[Hj and indefinite-metric linear spaces since the 1940's, [B]- In this context it corresponds to what is 
termed as "77-pseudo-Hermitian" in pQ, where 77 is an a priori fixed indefinite metric operator. The relevance 
of the indefinite-metric theories and PT-symmetric systems has been considered in [§]. The definition of 
a pseudo-Hermitian operator given in pp (and used below) is slightly different from the one used in earlier 
publications, e.g., [HI El 01]. As explained in detail in [TJ]|, this slight difference has important conceptual and 
technical ramifications. In particular, together with the idea of using biorthonormal systems pp it opens up 
the way for the construction of all possible metric operators, leads to the important observation that there 
is a positive-definite inner product rendering the Hamiltonian Hermitian for the cases that the spectrum is 
real [2J, and reveals the nature of the connection with antilinear symmetries such as PT, 0. 
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Hamiltonian operators that have a discrete nondegenerate spectrum. In particular, we will 
focus our attention mainly on the cases that the spectrum is not only discrete and nonde- 
generate but also real (and bounded from below). It is an operator with the latter properties 
that can serve as the Hamiltonian for a unitary quantum system, ^3] . If complex eigenvalues 
are present, we identify the vector space underlying the physical Hilbert space with the span 
of the eigenf unctions with real eigenvalues and restrict the Hamiltonian to this vector space 

dD mug. 

Let Ti be a given separable Hilbert space with inner product (-|-) and H : Ti — > Ti be a 
linear operator. Then H is called a pseudo-Hermitian operator pQ if there exists a Hermitian 
invertible operator 77 : Ti — > Ti satisfying 

H ] = riHr]- 1 , (1) 

where for any linear operator A : Ti — > Ti, stands for the 'adjoint of A\ i.e., the unique 
operator : Ti —>■ Ti satisfying (-\A-) = (A* ■ |-). The operator 77 entering the defining 
relation (fl|). which is sometimes referred to as a metric operator, is not unique ITT]. In 
fact the set Uh consisting of all metric operators is always an infinite set. A simple property 
of a pseudo-Hermitian operator is that it is Hermitian with respect to the possibly indefinite 
inner product (-, •)„ := (-\r)-), i.e., (-,H-) ri = (H-, •}„, [JJ. 

Next, suppose that H has a complete set of eigenvectors ip n G H, i.e., it is diagonalizable. 
Then one can construct the vectors <p n G Ti that together with ip n form a biorthonormal 
system for the Hilbert space, i.e., 

(0n|Vv) = S mn , ^ Mi&nl = 1- (2) 

n 

Using the properties of such biorthonormal systems, one can prove the following characteri- 
zation theorem [2]. 

Theorem: For a diagonalizable linear operator H with a discrete spectrum the fol- 
lowing conditions are equivalent. 

(cl) The spectrum of H is real. 

(c2) H is pseudo-Hermitian and the set Uh includes a positive-definite metric operator 
V+- 

(c3) H is Hermitian with respect to a positive-definite inner product (•>•)+, e.g., 
■= {'\V+-)- 

(c4) H may be mapped to a Hermitian operator h : TC — > Ti via a similarity transfor- 
mation, i.e., there is an invertible operator p :TL — > Ti such that 

h := pHp- 1 (3) 

is Hermitian. 
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If one (and therefore all of) these conditions hold, one has the following spectral resolutions 
for H and H'. 

H = Y J E n \lpn){<t>nl = J2 E n\<Pn}(lpn\. (4) 

n n 

Furthermore, a positive-definite metric operator r] + is given by 

V+ = ^2\<Pn)(<Pn\, (5) 
n 

and a canonical example of the invertible operator p whose existence is guaranteed by con- 
dition (c4) is p = y^+. 2 

The metric operator rj + plays the same role in pseudo-Hermitian quantum mechanics 
[TU] as the metric tensor does in general relativity, ^H]- It allows for the construction of the 
physical Hilbert space 7i p h y s and the observables of the system. The Hilbert space 7i p h y s has 
the same vector space structure as 7i but its inner product is given by 

(v)+ = (v)n + -=(-N+'>- (6) 

The observables O of the theory are linear Hermitian operators acting in 7i p h ys , P"3| - They 
can be obtained from the Hermitian operators o acting in 7i according to 

O = p- l op. (7) 

The formulation of the dynamics and the interpretation of the theory are identical with those 
of the conventional quantum mechanics. Pseudo-Hermitian quantum mechanics shares all 
the postulates of conventional quantum mechanics except that the inner product of the physical 
Hilbert space 7^ p h ys is not a priori fixed but determined by the eigenvalue problem for a linear 
(Hamiltonian) operator that acts on a reference Hilbert space 7i. 

As we mentioned above the formulation of the theory does not fix the reference Hilbert 
space H. For systems with a finite-dimensional state space, one usually identifies 7i with 
the complex Euclidean space, i.e., with usual Euclidean inner product: (i/j\4>) '■— ip* ■ <p, 
where a dot means ordinary dot product of vectors, jTT] . For PT-symmetric theories defined 
on the real axis, e.g., for H = p 2 + ix 3 , the natural choice for 7i is L 2 (IR), [TZj. However, 
for PT-symmetric theories that are defined on a complex contour T, such as H = p 2 — x 4 , a 
natural and useful choice for the reference Hilbert space 7i has not been available. The main 
purpose of this article is to offer a satisfactory resolution of this problem by showing how one 
can formulate and describe the same theories using equivalent PT-symmetric Hamiltonians 
whose eigenvalue problem is defined in L 2 (IR). This 'real description' facilitates the under- 
standing of the physical content of these theories. It allows us to use the usual mathematical 
tools of conventional quantum mechanics and deal with the manifestly non-Hermitian form 
2 For a mathematically rigorous discussion of pseudo-Hermitian operators, see |18| . 
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of the Hamiltonians such as H = p 2 — x A whose non-Hermiticity stems from their domain of 
definition rather than their explicit form. An alternative but less practical approach is to de- 
velop a pseudo-Hermitian description of PT-symmetric systems that is based on the choice: 
7i = L 2 (T). This 'complex description' clarifies the underlying mathematical structure of the 
formulation of PT-symmetric quantum mechanics that is based on the charge-conjugation 
operator H3 EEJ E3 • 

2 Moving Back to the Real Line 

Suppose JF is the set of real-analytic functions 3 ip : M —>■ C and H : T — > T is a linear 
operator of the form 

H =[p-A(x)} 2 + V(x), (8) 

where A, V : 1R — > C are piecewise real-analytic functions, pip(x) := —ii/j'(x) for all ip G J 7 , 
and a prime stands for a derivative. A particularly well-studied example is 

H = p 2 + x 2 (ixY, i/e(-2,oo). (9) 

The main observation that has led to the current interest in PT-symmetric quantum me- 
chanics is that for certain non-real choices of V (and A = 0), for example (jHJ) with v > 0, 
the operator H has a real and discrete spectrum provided that its eigenvalue problem is 
solved along an appropriate contour V in the complex plane, [4 . 4 This was a rather intrigu- 
ing observation because generically the operator H, which we will call the Hamiltonian, is 
manifestly non-Hermitian with respect to the L 2 -inner product. 

A typical physicist who is not familiar with the subject would immediately reject the 
statement that "if = p 2 — x 4 has a discrete spectrum." 5 Indeed, This statement is neither 
true nor false, because the eigenvalue problem for a linear operator defined on an infinite- 
dimensional vector space is well-posed only for specific choices of the domain of the operator. 
In the case of differential operators such as (jH)), in particular ©, the determination of 
the domain is related to the choice of the asymptotic boundary conditions. A nontrivial 
observation made in jl] is that one obtains a discrete spectrum for @ provided that one 
imposes the asymptotic boundary conditions along an appropriate contour T in the complex 
plane. 6 This means that one has to identify the eigenvalue equation for (jSJ) with its complex 

3 Notice that T n L 2 {R) is a dense subset of i 2 (R). 

4 A mathematically rigorous proof of this statement is given in \2'6\ I24| . 

5 This Hamiltonian corresponds to the choice v = 2 in ©. 

6 For the cases that v is an integer greater than —2, so that the potential term in @ is a monomial, 
this was known to mathematicians 26^. We will give a proof of the discreteness of the spectrum for all 
v G (— 2,oo) in the Appendix. 
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(holomorphic) extension 

r rl 1 2 

+ V(z)}V n (z) = E n * n (z), (10) 



and seek for solutions \I/ n such that 

|\l/„(z)| — > exponentially as z moves off to the infinity along T. (11) 

Note that the contour F is generally the graph of a (continuous piecewise) regular curve 
|77] parameterized by s G R, i.e., there is a (continuous piecewise) differentiable function 
C : R — > C with non-vanishing first derivative such that 

r = {c(*)|seR}, (12) 

and that lim^-toc, 9?[£(s)] = ±oo. Here and in what follows 3? and S respectively mean 'real' 
and 'imaginary part of. Clearly, we may state the boundary condition as 

l^n(C( s ))l — ¥ exponentially as s — > ±oo. (13) 

For the Hamiltonians 0, it is the choice of an appropriate contour T and the imposition 
of the boundary conditions (fTT]) that lead to a discrete set of nontrivial solutions for (fTDjh 
The same holds for various generalizations of ©, jSJ I2H]- I n general the contour T is not 
uniquely determined by the mathematical considerations, though it is required to stay in the 
so-called Stokes wedges in the asymptotic region, i.e., where s — > ±oo. In particular, there 
is a preferred choice for the asymptotic shape of T that maximizes the decay rate of the 
solutions of (JTUJ). This corresponds to the bisector of the appropriate Stokes wedge. Making 
this choice for the Hamiltonians we have [3] 

lim arg[C(s)] = -0± (14) 

where 'arg' abbreviates 'argument of and 

^ = -=2(^T4)' K'=*-Q*. (15) 

Next, we identify the real and imaginary axes of C with the x- and y-axes of the usual 
Cartesian coordinate system on M 2 = C, so that z = x + iy, and consider a general smooth 
contour T such that §t[T(x+iy)] is an increasing function of x := 9ft(z). 7 Then we can express 
the function ( in terms of a differentiable real-valued function / : R — > R according to 

C(x)=x + if(x). (16) 



7 This is not a strong condition. One can always choose such a contour for the purpose of defining 
boundary conditions 
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The condition that ( is a regular curve is also satisfied, because |£'(x)| 2 = 1 + f'(x) 2 ^ 0. 

Now, we wish to restrict the complex differential equation (fTTH) to the contour T, and 
obtain an equivalent real differential equation with generally complex coefficients. Along V 
we have z = ((x) = x + if(x). A simple change of variable z — > x + if(x) in (fTUj) yields 



-g(x) 2 



ia(x) 

ax 



" • i<iuY\r"(.r) r 



ia(x) 

ax 



+ v(x) > ip n {x) = E n ip n (x), (17) 



where 



g(x) := [ax)]- l = [l + if{x)}-\ a{x):=g{xY 1 A[x + if{x)l (18) 
v{x) := V[x + if{x)\, ii n {x):=^ n [x + if{x)}. (19) 

The complex differential equation (fTTH) together with the boundary condition (fTT|) (alterna- 
tively (|T3|) ) is clearly equivalent to real differential equation (fT7|) together with the boundary 
condition 

|?/>(x)| — > exponentially as \x\ — > oo. (20) 

The analyticity properties [23] of and consequently of ip n together with the condition ([2*U|) 
implies that ip n G L 2 (M). In other words, the eigenvalue problem for the Hamiltonian (jH|) 
defined by Eq. (jl(Jj) is equivalent to the eigenvalue problem for the Hamiltonian 

H' := g(x) 2 [p - a(x)} 2 - g(x) 3 f"(x)[p - a(x)\ + v(x). (21) 

viewed as an operator acting in L 2 ( 



3 Consequences of Imposing PT-symmetry 

Let £ : K. — >• C be a function. Then under the joint action of the parity P and time-reversal 
T operators, C,(x) — > PT £(x) PT = £(— x)*. Applying this rule to the Hamiltonian (|21j) and 
using PT p PT = p, we find 

PTH'PT = g{-x)* 2 \p-a(-x)*] 2 -g(-x)* 3 f' , (-x)*\p-a(-x)*]+v(-x)*. (22) 
In particular, demanding H' to be PT-symmetric yields 

g(- x y 2 = g(x) 2 , g(- x yf»(- x y = g (- x y f »( x ), (23 ) 

o(— x)* = a(x), v(—x)* = v(x). (24) 

In view of Eqs. fITHj) . (fTDj) . (|2^j) . and (j2U), the fact that / is a real- valued function, and x 
takes zero as a value, we have 

/(*) = /(-x), A ( u T\ u =-[x+if(x)] = A[x + if(x)}, V(uy\ u= _ [x+mx)] = V[x + if(x)}. 

(25) 
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The first of these equations imply that along the contour T, z(—x)* = —z(x). Therefore 
the condition that H' be PT-symmetric implies that T has reflection-symmetry about the 
y- (or imaginary-) axis. The second and third equations in (J23j) and the assumption that A 
and V may be analytically continued onto the contour V indicate that they are separately 
PT-symmetric, i.e., 

PTA(x)PT = A(x) } PTV(x)PT = V(x). (26) 

These are equivalent to requirement that the original Hamiltonian (jSJ) be PT-symmetric. 

In summary, the Hamiltonian (jHJ) and the contour T are PT-symmetric if and only if the 
Hamiltonian ()21|) is PT-symmetric. In the following we will only consider the cases that 
these conditions hold. 



4 Wedge-Shaped Contours 

The simplest possible PT-symmetric choices for the contour T are the wedge-shaped con- 
tours: 

T(x) = x[l — i sign(x) tan 9], (27) 

where sign(x) := x/\x\ for x ^ 0, sign(0) := 0, and 9 G [0,7r/2). Clearly V is not a regular 
curve at x = 0. Therefore, we will smoothen it in a small neighborhood of x = 0, say 
according to T — > T e , where 

r e (a;) := x + if e (x), (28) 

, — \x\ tan# for \x\ > e . . 

fe(x) ■= < ' . T (29) 



(fe(x) for |x| < e, 



etan6* 



ip e [x) 



x\ „ fx 



, Hi. 2 



(30) 



and e G M + is an arbitrary constant. Note that f e is a twice-differentiable function that can 
be substituted for / in the expression (|21|) for the Hamiltonian H' and that its maximum 
value is / e (0) = — 3etan#/8. Figure [T] shows a plot of f € . Furthermore, in view of and 
P?|l . we have 

f sec(9)e- idsi ^x for Id > e 

r e (x) = < (3i) 

I x + i(p t (x) for \x\ < e. 
In what follows we shall consider the contours of the form ()31|) which yield the wedge- 
shaped contours (}2"Tj) in the limit e — >• 0. 

Setting / = f e in (fTHj) and using (j2T?j) . we obtain 

for |ar| > e : f'(x) = -tan(0) sign(x), f"(x) = 0, g(x) = cos{9) e iesign{x \ (32) 
for|x|<e: f'(x)=<p' e {x), f"(x) = <p»(x), g(x) = le (x), (33) 
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Figure 1: Plot of y = / e (x). f e has a maximum at x = with 
value / e (0) = — 3e tan 9/8. The angle 9 is also displayed 



where 



76 (z) 



tan# 
2 

3 tan# 



2e 



!) 3 



i + 



itan# 



These relations together with (jI8}l. (fT9"j) . and (|2~TJl then yield 



#' = i/i e) +Ht + H 



(<0 



where 



A 



dar|a;)(ar|, A_ := / dar|a;)(a;|, A e := / dx\x)(x\, (38) 

#± := cos 2 (0)e ±2ie {p-sec(^)e Tie A[sec(^)e^x]} 2 + F[sec(0)e Ti0 x], (39) 
# e := le (x) 2 [p - a e (x)} 2 - le {xfip" e {x)\p-a e {x)\ +v e (x), (40) 
a e (x) := 7 e (x) _1 A[x + i<^ e (rr)], v e (x) :—V[x + i(p e (x)]. (41) 

Note that PTA^PT = A W and PTA e PT = A e . These together with and (IHUl) - 
(|41|1 yield the following relations that are clearly consistent with the PT-symmetry of H'. 



(0 .. 



(6) 



(34) 
(35) 
(36) 



:= A E H e A ej (37) 



PTPi e) PT = H [ j\ 



PT H f PT = H f . 



(42) 



In practice, to solve the eigenvalue problem for H', we may solve the corresponding 
differential equation for |x| > e in the limit e — > and match the solution at x = by 
enforcing appropriate continuity requirements. As we shall see below the latter yield a pair 
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of boundary conditions at x = 0. It is the Hamiltonians H± together with these boundary 
conditions at x = and the requirement: ip n G L 2 (R) that determine the eigenvalues E n . 

The Hamiltonians H± take a simpler form in terms of the scaled position and momentum 
operators: 

x 

x := -, p:=cos#p. (43) 

cosU 

The classical analog of x corresponds to the arc- length parametrization of the contour T, 
[27]. Using 1|BS |) and we have 

H± := e ±2ld [p - e^ 9 A(e^x)] 2 + V (e Tie x). (44) 

The boundary conditions at x = may be obtained by integrating both sides of the 
eigenvalue equation for H' over the interval [— e, e] and taking the limit e — > in the resulting 
expression. Doing an integration by parts, using the fact that A and V are continuous 
functions, noting that 



Ve (±e) = ^'(±e) = 0, <^(±e) = 1, 7e (±e) = (1 =F itani 
and introducing the notation 

^n(0 ± ) := lim .tl> n (x), ^(0 ± ) := lim ^(a;), 

x^O* X— >0 ± 

we find the following boundary condition at x — 0. 

Imposing the condition that ip n be continuous at x = 0, i.e., 

^(0 ± )=^(0), (46) 

reduces (|43)) to 

e- 2 ^(0-)=e 24 V;(0 + ), (47) 

or equivalently to 

K(0-)| = K(0 + )| and (48) 
arg[<(0-)] = arg[<(0-)] + AO if <(0±) ^ 0. (49) 

Therefore, for ip n to be different iable at x = either = or i/j' n (0) = 0. 
For a PT-invariant eigenf unction ip n , where 

Ms) = Mx)*, = (50) 
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and in particular 



^n(O-) = ^n(0 + )=^„(0)GM (51) 

<(0") = -<(0 + )*, (52) 



(gnj) implies that 



either <// (0") = ^(0 + ) = or arg[^(0 ± )] = | =f 20. (53) 

As a result is differentiable at x = if at least one of the following conditions hold: 
1. t/>'(0) = 0; 2. 9 = and ^'(0) is imaginary. 8 

Having derived the explicit expression for the boundary conditions at x = we can 
identify the eigenvalue problem for the initial Hamiltonian H and the contour ([27]) with that 
of 

H' = Ai 0) tf_Ai 0) + AfiJ+Af (54) 

and the requirement that the eigenfunctions belong to L 2 (1R) and satisfy the boundary 
conditions (jjB]) and ([47)1 . For real eigenvalues, where we may choose to work with the 
PT- invariant eigenfunctions, we have the boundary conditions (|51|). ()48p. and ()53jl. 

5 Application to H = p 2 + x 2 (ix) u 

For the Hamiltonians we have 

A(x) = 0, V{x)=i v x» + \ e = e v ,= -^—. (55) 

Inserting these in (|44|) . we are led to the following remarkable result. 

#± = e ±2 ^ H v+ 2, (56) 

where 

#jv : = p 2 + jxl^ for JVGK. (57) 
Therefore, in view of (|54j) . we have 

H' = e - 2ie »A {0) H u+2 Ai 0) + e 2i ^Af H v+2 Af . (58) 

The eigenvalue problem for the Hamiltonian © that is defined by the contour (|2*7jl with 9 
given by (jB3j) is equivalent to the eigenvalue equation 

e 2^si g n(x) [_^( x ) + | x ^+ 2 ^ n ( x )] = Eifj n {x) for x ^ 0, (59) 

8 In conventional quantum mechanics, where = 0, the PT-symmetric eigenfunctions of a PT-symmetric 
Hamiltonian of the standard form p 2 + V(x) are either real and even (where condition 1 holds) or imaginary 
and odd (where condition 2 holds). For an example see 1 1 3 j . 
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where if) n are required to be continuous elements of L 2 (R) satisfying 

e-*V„(0-) = e*V B (0 + ). (60) 

Next, we show that the eigenf unctions tp n never vanish at x = and they are necessarily 
non-differentiable at this point. 9 

Lemma: Let ip n G L 2 (R) be a continuous solution of (JHHJ) and (jHUj) with z/ > —2 and 
7^ and ^> n ± : K U {0} — > C be its restrictions: ?/vt(x) := ^n(x) for all ±x G M + , 
V> n± (0) := Vn(0), and < ± (0) := <(0±). Then 

^/O/^). (61) 

Proof: Clearly (J5T?|) and (JBTlj) are respectively equivalent to 

-<±(x) + |x|^Vn±(x) = e^" £ n ^n±(x) for ± x G R ± , (62) 
e- 2 ^<_(0) = e 2 ^< + (0). (63) 

Multiplying both sides of (J62j) by ipn±i integrating over M ± U {0}, and performing an 
integration by parts yield 

±Vn±(0)*V/ ± (0)+ || ^ ± HI + || W^Vni 111= || Vn± 111, (64) 

where for all f± : M ± -> C, || f ± |||:= / R± |£ ± (x)| 2 dx. 10 Now, if at least one of ^(0), 
ijj'(0 + ), and ip'(0~) vanishes, then so is the first term in (J64)) . This implies that e T "E n 
must be real for both choices of the sign. For v > —2 and v ^ 0, this is only possible 
if E n = 0. But then the right-hand side of ()64|) vanishes, while its left-hand side is 
strictly positive. This is a contradiction proving (|6ip. □ 

A direct implication of ()61|) is that if v > —2 and v ^ 0, then for all n, ip n fails to be 
differentiable at x = and that we can always normalize ip n so that ip n (0) = 1. 

For real eigenvalues E n we can take ip n to be PT-invariant and for the cases of interest, 
namely v > 0, the boundary conditions on the eigenvalue equation (|59|) take the form 

^n(O-) = ^n(0 + ) G E, (65) 

K(0")| = K(0 + )|, (66) 



9 V>n is necessarily twice differentiable at all x ^ 0. 

10 Note that because |x| 2+l/ is bounded from below, || ip' n± |||. and || |x|"/ 2+1, n ± ||j_ are finite numbers, 
§10.1]. 
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An interesting particular example is the Hamiltonian 

H = p 2 - x 4 , 



(68) 



which corresponds to u = 2 and 

H ± = e ±! i [p 2 +x 4 ], (69) 

with eigenvalue equation 

ef sign(x) j_^ (x) + x 4^ (x) ] = for x _^ 0? (7Q) 

and boundary conditions (JUSJ), (pjj) and 

axg^^)] = (71) 

The switching of the sign of the potential term from minus in ()68|) to plus in ()69|) and 
(J7UJ) is quite remarkable. As seen from ([57)1. (|58j) and ([59*]). this is a characteristic feature 
of the Hamiltonians H of the form Q. In view of the discreteness of the spectrum of the 
Hamiltonians for N > 0, j2H], this phenomenon provides invaluable insight in the origin 
of the discreteness of the spectrum of H. Indeed, as we shall show below, it leads to a 
rigorous proof of the fact that for all v 6 (—2, oo) the spectrum of H is discrete. Note 
that here and in what follows the spectra of H^, H', and H are respectively defined by the 
exponentially vanishing boundary condition at ±oo along M, the latter together with the 
boundary conditions (|63|) - (}6T|) at x = 0, and exponentially vanishing boundary condition 
at ±oo along the contour (}3*T]) with 9 = 9 U . 

To establish the discreteness of the spectrum of H' (and consequently H), we use the 
equivalence of the eigenvalue problem for H' with Eqs. (J62|) and (J63)) . and note that in terms 
of the functions y± : [0, oo) — > C defined by 

y±(x) :=Vn±(±x), (72) 

(p)2*j) takes the form 

-j/£(x) + x u+2 y±(x) = A±3/±(x), for x G [0, oo), (73) 

where 

A ± = e^ 2 ^E n = e T S£„. (74) 

The eigenvalue problem for H' is equivalent to finding the solutions y± of (J73|) that belong 
to L 2 [0, oo) and satisfy 

y_(0) = y+(0)^0, (75) 
2/-(0) = -e^ + (0)^0. (76) 
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This problem may be treated using the classical theory of singular boundary- value problems 
developed mainly by Weyl, [23 §10]. In the Appendix, we will use some basic results of this 
theory to give a proof of the discreteness of spectrum of H for all v G (—2, oo). 

We close this section by pointing out that the formulation of the eigenvalue problem 
for H as the differential equations (fTHj) with boundary conditions fjTHJ) and (fTBJ) is also of 
practical importance because it allows for the immediate application of the known numerical, 
perturbative, and variational methods that are tailored to deal with functions of a real 
variable, [HO]- It should also be interesting to see if one can obtain an alternative proof of 
the reality of the spectrum using this formulation. 



6 Square Well Placed on a Wedge-Shaped Contour 

Consider the Hamiltonian H = p 2 + V(x) for the ordinary Hermitian infinite square well 
potential 



V(x):=\ (77) 
[ oo tor \x\ > -j, 

where L G IR + . If one solves the eigenvalue problem for this Hamiltonian on the real axis 
one finds an infinite discrete set of eigenvalues 

E$» = n G Z+. (78) 

As this Hamiltonian is both Hermitian and PT-symmetric, one may choose to work with 
normalized PT-invariant eigenfunctions which are given, up to an arbitrary sign, by [T3] 



X 1 

Tin I — I — 
L 2 



1 + t-iY 

»n := \ ' ■ (79) 



We wish to explore the consequences of defining the eigenvalue problem for the square well 
Hamiltonian using a wedge-shaped contour (|27|) with arbitrary angle 9 G (0,7r/2). n 

Pursuing the approach of Sec. 4, we find that the eigenvalue problem for this system is 
equivalent to the following boundary-value problem. 

-< ± (x) = e^ e E n ^ n± (x) for ±xG[0,^], (80) 

^ n _(0) = VW(0), e- 2j V;_(0)=e +2 ^ + (0), (81) 
^±(±|) = 0. (82) 



11 Taking the v — > oo limit of JSJ) one obtains a similar square well Hamiltonian (with L = 2 and = 6 V — > 
7r/2), For large but finite value of v this Hamiltonian has an infinite number of positive real eigenvalues 
all of which are proportional to v 2 . Therefore in the limit v — > oo, real part of the spectrum is mapped to 
(the point at) infinity. The spectral problem considered in this section is different from the one treated in 
|31j. for we view the potential l|77|l as given and take as a free parameter. We will see that for large 
{9 > 7r/4) the spectrum is entirely complex. 
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Clearly ifj n ± determine the eigenfunctions ip n of the system according to 

V>„(x)=V„±(x), if ±xe[0,f]. (83) 

They belong to 

H' := {V e L 2 [-§, f ] I V<±3r) = 0} . (84) 

The eigenvalue problem (jHUj) - (JH2j) may be easily solved: Zero is an acceptable eigenvalue 
only for 9 = 7r/4. The corresponding PT- invariant eigenfunction is given by 

V>(x) = ±c (x =f | ) for ±xe[0,f], (85) 

where c is a real normalization constant. The eigenfunctions with nonzero eigenvalues have 
the form 

V„±(x) = c± e^ ±x + d ± e"^ ±x , (86) 

where u n ± := e Tie \/E n ~ and 

c- = \ [(1 + e 2 ^)c + + (1 - e™)d + ] , d_ = i [(1 - e 2 *V + (1 + e 2 *V+] , (87) 

Eqs. (|57|) and (jHHj) follow from the boundary conditions (|ST|) and (jH^j) . respectively. They 
have a nontrivial solution provided that the eigenvalues E n satisfy a transcendental equation 
that takes the following simple form in terms of the variable u n := cos(8)Ly/E^, 

tan(#) sinh[tan(6 l )-u n ] = sin(w n ). (89) 

For 6 = it reduces to sin(w n ) = 0, and one recovers E n = En . But for 6 > it has a 
finite number N(9) of real solutions where iV is a decreasing function of 9. In particular, for 
6 > 7r/4, N(9) = and there is no real solution. As one decreases the value of 9 from 7r/2 
down to zero one encounters an infinite strictly increasing sequence {Se} of exceptional points 
|32j . The angles 9 for the corresponding wedge-shaped contours form a strictly decreasing 
sequence {9i} that converges to zero. Table [U lists the values of the first five exceptional 
points and the corresponding angles 9g_. 

In general, the number of real eigenvalues are given by 

{2t-\ for 9 e (9 t +i, 9 1) with i > 1 
2£ - 2 for 9 = 9 e with i > 2 (90) 
1 for 9 = 9 1 = it /A. 

Because the eigenvalues are nondegenerate, the dimension of the invariant subspace spanned 
by the eigenfunctions with a real eigenvalue is N{9). This A^(6 l )-dimensional subspace is the 
underlying vector space V for both the reference Hilbert space (7i) and the physical Hilbert 
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i 


1 


2 


3 


4 


5 







61.58 L~ 2 


200.9 L~ 2 


418.9 L~ 2 


715.7 L~ 2 


e e 


45.00° 


14.81° 


9.88° 


7.59° 


6.23° 



Table 1: The first five exceptional points Eg and the 
corresponding exceptional values 9g of 9. 

space (7i p hys) of the system. For 9 = 0, N(9) = oo and TL, 7i p h y s, and W coincide. But 
for 9 > 0, V is finite-dimensional. In particular, for 9 > 9\ = tt/4 the vector space V is 
zero-dimensional and the system does not admit a unitary quantum description. 

Another peculiar feature of this system is that the dimension of the physical Hilbert space 
takes even values only for the exceptional values 9g of 9 with I > 2. As these constitute a 
measure zero subset of [0, 7r/4), the physical Hilbert space is generically odd-dimensional! 

Three comments are in order. 

1. If one defines the eigenvalue problem using the Neumann boundary conditions at x = 
±L/2, i.e., requires tp' n ±(±L/2) = 0, the (nonzero) eigenvalues are given by Eq. (}89|) 
with the sign of the right-hand side changed. The corresponding pseudo-Hermitian 
quantum system shares the general features of the square well system discussed above. 
The only difference is that for all values of 9, zero is an eigenvalue with a constant 
eigenf unction. In particular the physical Hilbert space is finite-dimensional for < 
9 < vr/2, infinite-dimensional for 9 = 0, and one- dimensional for it /A < 9 < ir/2. 

2. The quantum system corresponding to the square well Hamiltonian placed on a wedge- 
shaped contour defines a PT-symmetric quantum system which is fundamentally dif- 
ferent from the PT-symmetric square well studied in |3*3*1 |33J • The latter system 
involves a non-Hermiticity parameter Z G [0,oo). As one increases the value of Z 
(starting from zero) one encounters an infinite sequence of exceptional points which 
correspond to a strictly increasing sequence {Zi\ of exceptional values of Z. As a 
result unlike the system introduced above the physical Hilbert space is always infinite- 
dimensional. In particular, for < Z < Z x the reference Hilbert space 7i coincides 
with W. 

3. For the square well system defined on a wedge-shaped contour, 9 = 0-- which 
corresponds to the Hermitian limit of the problem - - is an accumulation point of 
the exceptional values 0£ of 9. This is the reason why for all positive values of 9 the 
physical Hilbert space is finite-dimensional. 12 This observation shows that changing 

12 It is not difficult to see that the same holds for negative 6. 
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the domain of the definition of a Hamiltonian from the real line to a complex contour 
can lead to completely different quantum systems. For example, for 9 = 9 2 the physical 
Hilbert space is two-dimensional. Therefore it describes the interaction of a spin-half 
particle with a magnetic field [33]. In contrast, for 6 = 0, the system describes the 
one-dimensional motion of a particle that is trapped between two impenetrable walls. 



7 Application of Pseudo-Hermitian QM for Square Well 
along the Wedge-Shaped Contour with 6 = 62 

The largest value of the angle 9 that corresponds to a nontrivial unitary quantum system 
is 9 = 62 ~ 14.81°. For this choice of 9 the Hamiltonian has two real eigenvalues. They 
are E\ ~ 9.09L -2 and Ei = £2 ~ 61.6L -2 . The corresponding eigenfunctions tpi and fa are 
given by (jHHj) and (jEHJ) where 

c n ± = c±(E n ), d n ± = d±(E n ), (91) 

Q±(£) := e^LVE, (93) 

and jv(i£) is an arbitrary real normalization constant. Substituting ()91|) and (}9~2"j) in (J8fi)l and 
using (JH3J), we have 

M*) - **M = ^ Si " [ "^"' 3 — f)] for ±x 6 [0,f], (94) 

sin ^ nj 

where Af n G IR + are normalization constants, and n = 1,2. 

The underlying vector space V for the reference and the physical Hilbert spaces is the 
two-dimensional subspace of TC' spanned by fa and fa. The reference Hilbert space TC is 
obtained by endowing V with the subspace inner product (-|-) induced from TC'. Choosing 
the normalization constants as a/"i ~ 1.226L _1 / 2 k and N2 ~ 0.717L -1 ' 2 ft, for some k G K + , 
we have 

<^i|-0i> = <^|V>2> = k\ {fa\fa) = (fa\fa) = tk\ (95) 

where 

r « 0.068. (96) 

Clearly {fa, i/j 2 } form a non-orthogonal basis of TC. We can use the Gram-Schmidt procedure 
[3E] to construct an orthonormal basis {£i,£ 2 } according to 

£1 := re Vi, e 2 := — ===. (97) 

KV 1 — r z 
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In this basis the Hamiltonian is represented by the following manifestly non-Hermitian 2x2 
matrix 

\ E 2 J \ 61.6 J 

We can compute the adjoint of H using its matrix representation (J98)) and determine its 
eigenvectors cj) n that together with ip n form a biortho normal system for TC. This yields 

Now, we are in a position to compute the metric operator rj + . In view of (J5J) and (|99J1 . it has 
the following matrix representation in the orthonormal basis {61,62}- 

1 -0.068 



0.068 1.009 



VT^F 1 l-r 2 

In view of this relation, we have, for all £, £ G V, 

r(£C2 + £d) , (l+r 2 )^C 2 



(e,c>+:=<eh+c> = « 



-2 



SIC: 



VT^T* 1 -r 2 

~ ^ 2 KCi- 0.068 (£C2 + £G) + 1-009 e 2 *C 2 ], (101) 

where £ n = (e n |£)> Cn = ( £ n\()i an d w = 1, 2. Note that the coefficient k~ 2 is a trivial scaling 
of the inner product. 

If we use ()101|) to compute the inner product of the eigenvectors i[> n , we find that as 
expected {tp%, ^2} form an orthonormal basis of the physical Hilbert space, (ip n , ip m )+ = 8mn 
for m,n = 1,2. This also shows that the Hamiltonian viewed as acting in 7i p h y s is a Hermitian 
operator. 

Next, we construct the physical observables O of the system. This requires the com- 
putation of p = y/rj + - The matrix representation of p in the basis {61,62} has the form 

,— . / 0.999 -0.034 \ , s 

p=^l^K- l [ . 102 

\ -0.034 1.004 J 

According to Q, the physical observables are given by O = ^£=0^^ wnere are 
arbitrary constants, S is the identity operator acting in TC, for £ = 1,2,3, are defined 
through their matrix representations in the basis {61,62} according to 

S £ = p^aip, (103) 



and <j£ are Pauli matrices. Specifically, 

^0 = : = , hi 







-c 


:)• 


0.068 


-1.007 


0.998 


-0.068 



1.005 
0.995 

1.002 -0.068 
0.068 -1.002 
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Using these relations and (J98j). we can show that indeed 

H » L~ 2 (35.5 S + 1.78 £ x - 26.2 E 3 ) . (104) 

Next, we compute the Hermitian Hamiltonian h of (|3*|) that is associated with i/. We 
can obtain the matrix representation h of h in the basis {£1,62} using either of (jlU2j) and 
flSEJ or ([TMil and ©. Both yield 

~ „ / 9.15 1.78 \ „ . . . 

h^L- 2 \ = L- 2 35.5 a + 1.78 cri - 26.2 a 3 ). 105 

\ 1.78 61.5 J 

Therefore, 

/i « L- 2 (9.15 |ei)(ei| + 1.78 (|ei)(e 2 | + |e a )<ei|) + 61.5 |e 2 )<e a |) . 

Having obtained the biorthonormal system {\ip n ), \<f>n)}, we can also compute the gener- 
alized parity V, time-reversal T, and charge-conjugation C operators of ^7j, namely 13 

(106) 

I, (107) 

(108) 

where * is the complex-conjugation defined by 

2 2 

* |C) := £>«IO*k«> = £<Ckn> kn), for all C e W. (109) 

n=l ra=l 

In particular, in the basis {61,62}, * is represented by ordinary complex conjugation V of 
complex vectors, 

*z:=z>, where z—\ ' | £ C 2 (110) 

V (e a |C> ) 

As explained in ^7j, unlike C which is always an involution (C 2 — 1), V and T need not 
be involutions. Requiring them to be involutions restricts the choice of the biorthonormal 
system. In the case at hand, this restriction amounts to fixing the normalization constant 
for the eigenvectors ipn as 

«=(l-r 2 ) -1/4 w 1.001. (Ill) 

Making this choice, we find that the matrix representations of V, T , and C, in the basis 
{61,62}, are respectively given by 



V 


■= \<f>i)(<h\ 


- \<f>2)(<fo\: 


r 


■= |0l}*(<? 


h\ ~ |02) * 


c 




- 1^2) (<h\ 



, VT^r 2 -r \ ( 0.998 -0.068 . 

V = ; « , 112 

-VT^ J \ -0.068 -0.998 

, 1 —£=s \ ( 1 -0.136 , 
T = P,. C= ( q ^ U [ _ i ]• (113) 
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See also [S3- 
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Using these relations we can directly check that indeed 

V 2 = T 2 = C 2 = 1, C = r^P, [H, C] = [H, VT] = 0. (114) 

In view of the identity VT = the PT-symmetry of H corresponds to the fact that H is 
a real operator with respect to the complex-conjugation ()109j) . i.e., * H * = H . An explicit 
manifestation of the latter relation is that if is a real matrix. 14 

8 Formulation Based on the CVT-Inner Product, Dis- 
cussion, and Conclusion 

In this paper we have presented a formulation of PT-symmetric theories defined along a 
complex contour in which the state vectors belong to the familiar Hilbert space of square- 
integrable functions. This formulation has a number of advantages. Firstly, it yields the 
necessary means for a straightforward application of the results of the theory of pseudo- 
Hermitian operators. Secondly, it provides a novel description of the Hamiltonians of the 
form (jUJ) that reveals the origin of the discreteness of their spectrum. Finally, it is practi- 
cally appealing for it allows for a direct application of the standard approximation schemes 
developed for solving differential equations on the real line |3Uj . 

In order to elucidate the practical aspects of our method we have considered the PT- 
symmetric system obtained by placing an infinite square well potential on a wedge-shaped 
contour T. We have conducted a comprehensive study of this model showing that as soon as 
one makes the characteristic angle 9 of the contour V different from zero (i.e., moves off the 
real axis) the physical Hilbert space of the system becomes finite-dimensional. The dimension 
of this space depends on 9. It changes at certain critical values of 9 that correspond to the 
exceptional spectral points associated with the system. The simplest nontrivial case occurs 
at the second exceptional point where 9 ~ 14.81° and the physical Hilbert space is two- 
dimensional. For this case we showed how one could employ the constructions developed in 
the framework of pseudo-Hermitian quantum mechanics to determine the explicit form of the 
inner product of the physical Hilbert space, the physical observables, and the corresponding 
Hermitian Hamiltonian. 

The results reported in this paper show that PT-symmetric quantum mechanics is in- 
deed a special case of pseudo-Hermitian quantum mechanics. In order to apply the pseudo- 
Hermitian quantum mechanics to PT-symmetric systems defined on a complex contour, 
one may employ the fact that these systems admit a convenient description in terms of 

14 Because the matrix representation H of the Hamiltonian is not symmetric, the definition of observables 
proposed in [SSj cannot be employed, |55|. 
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PT-symmetric Hamiltonians defined on the real line. The latter may be treated most per- 
spicuously within the framework of pseudo-Hermitian quantum mechanics. In particular, 
one can compute the observables of the theory and explore its classical limit as outlined in 

mum. 

There is also a more direct, but less practical, pseudo-Hermitian description of PT- 
symmetric systems defined on a complex contour T. This is also suggested by the analysis of 
Sec. 2. 15 It involves identifying the reference Hilbert space Ti with L 2 (T), where the contour 
T is viewed as a one-dimensional real submanifold of IR 2 = C, i.e., a continuous (piecewise 
regular) plane curve. The relationship between this 'complex pseudo-Hermitian description' 
and the 'real pseudo-Hermitian description' that is based on transforming the system onto 
the real line may be reduced to the action of a diffeomorphism Q of the complex plane that 
maps the real axis onto the contour T. This mapping may be identified with the arc-length 
parametrization of T. In view of we can parameterize T by the ^-coordinate. We can 
use this parametrization to define the arc-length parameter: x = T{x) := f* a/1 + f'(s) 2 ds. 
Note that for the contours of interest JF : R — > E is a diffeomorphism. The restriction of Q 
onto the real axis defines the following mapping of IR onto T. 

g(x):=x + if(x)=F- 1 (x)+if(F- 1 (x)), for all x G K. (115) 

This in turn induces a unitary operator u g : L 2 (R) — > L 2 (T) defined by 16 

(u g tfj)(z) :=i){<g-\z)), for all ip G L 2 (U), zeT. (116) 

Alternatively, setting ^ := u g ij) we have 

$(z) = ^(x) if and only if z = £(x). (117) 

The statement that u g is a unitary operator means that for all ip, <fr G L 2 (IR) 

(^M)r=(V#>, (118) 
where (-|-)r is the inner product of £ 2 (r), i.e., 

(tt|$) r := J ^(zy$(z)dz. (119) 

The validity of Eq. ()118|) becomes obvious once we identify T with a plane curve and view 
the right-hand side of ()119|) as a line integral. Letting \I/ := u g tp and $ := u g <p and using 
(fTTBT) and (fTToT). we have 

{u g i)\u g <P) T = (tf|$) r = / ^(£(x))*$(£(x))dx= f ^(x)>(x)rfx= 

Jr Jr 

15 See also [231 

16 One might try to express u g in the form e 1 ^ ^' 9 ^ 2 for some complex-valued function G by extending 
the results of |3H] . 
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An important property of u g is that it establishes a one-to-one correspondence between 
the ingredients of the two pseudo-Hermitian descriptions of the system; to each linear op- 
erator A acting in L 2 (R) it associated a linear operator A := u g Au~ l acting in L 2 (T). In 
particular, it maps the charge-conjugation operator C := r]~ x P of the real description to the 
charge-conjugation operator C : L 2 (T) — > L 2 (T) of the complex description according to 

C:=u g Cu-\ (120) 

In view of the results of for the Hamiltonians (jOJ with v > 0, the operator C is nothing 
but the charge-conjugation operator introduced in [20] ■ In fact, what the authors of [20] do 
is to define C on the real line (though they use the same symbol for both C and C) , perform 
the diffeomorphism u g to obtain C, and then use it in a contour integral along T to define 
their CVT- inner product: 

(*,$)cpt:= j[CVT^{z)]<S>{z)dz for f,$6i 2 (r). (121) 
Note that in the real description [T7j . 

(V, 0)cpt := / [CPT^(x)]0(x) C /x = (i/ilvM = (il>, 0> + for if,, G L 2 (R). (122) 



Moreover, the eigenf unctions ^ (respectively ip n ) 17 form an orthonormal set with respect 
to (-,-)cvt (respectively (-,-)+), 

($ m , ^ n )cvr = 5 mn = (ip m , ip n ) + = (ip m \r] + \ip n }. (123) 

Next, we introduce a metric operator rfl : L 2 (T) — > L 2 (T) and the corresponding inner 
product (-, : L 2 (T) x L 2 (r) -> C according to 

i£ := ti^ 1 , (■,->S:=(-|tf->r. (124) 

In view of the identity \I/„ = u g ip n and the fact that u g is unitary, we then find 

Equations (|123|) and (|125|) show that \l/ n , which are supposed to form a complete set, are 
orthonormal with respect to both the CPT-inner product (|121|) and the inner product (■,■)+. 
This proves that these two inner products are identical. Therefore, the formulation of PT- 
symmetric quantum mechanics based on the CPT-inner product, as outlined in [20], admits 
a complete description in terms of the theory of pseudo-Hermitian operators. 

17 Recall that according to the analysis of Sec. 5, the eigenfunctions ip n and ^> n are related via \P n (C7(x)) = 

-0n(x). 
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Note: After the completion of this project, I discovered a preprint of Znojil |H] where 
he considers the analytic continuation of the PT-symmetric square well of Ref. [33] onto 
a smooth complex contour. The spectral properties of this system is similar to the one 
considered in Sec. 6. In both cases the spectrum is determined through a set of boundary 
conditions at the intersection point of the contour and the imaginary axis. The main dif- 
ference between the two systems is that the defining boundary conditions used in [JT] are 
postulated whereas those used in Sec. 6 are derived. As explained in Sec. 4, the latter are 
the general boundary conditions associated with the wedge-shaped contours. 
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Appendix: Discreteness of the Spectrum of (E3) 

Theorem: The spectrum of the Hamiltonians H = p 2 + x 2 (ix) u defined by the contour 
fj27jl with 9 = 6 U := ttu/[2(u + 4)] is discrete for all v G (-2, oo). 

Proof: For v = this statement is well-known to hold j2H]- To prove it for v ^ 0, we 
prove the equivalent statement that for all v G (—2, oo) the following boundary-value 
problem has a solution only for a discrete set of values of E n . 18 

-yi(x) + x» +2 y ± (x) = X ± y ± (x) for x G [0, oo), (126) 

A ± = e T2W »E n G C, (127) 

y ± GL 2 [0,oo), (128) 

j/_(0) = y+(0) ^ 0, (129) 

^( ) = -e 4 *SV(0)^0. (130) 

Let A G C be arbitrary, and consider finding solutions y(-; A) of 

-y"(x) +x u+2 y(x) = \y{x), with v > -2, xg[0,oo), (131) 

that belong to L 2 [0, oo). Then because x v+2 is bounded below by zero, one has the 
so-called limit point case [23 §10.1] where there is at most one linearly independent 
L 2 -solution and such a solution exists for all non-real A and has the form 

y(x; A) = C(\)[ yi (x; A) + m(A)y 2 (x; A)], (132) 



18 The equivalence of this statement with that of the above theorem is established in Sec. 5. E n are the 
eigenvalues of H. 
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where C(A) G C — {0} is a constant, y\ and ?/2 are the fundamental solutions of (|13H) 
satisfying 

yi(0;A)=0, 2/i(0;A) = -l, y 2 (0;A) = l, ^(0;A)=0, (133) 
and m : C — > C is a function having the property [23 §10.2] 

m(A*) = m(A)\ (134) 

Now, consider the boundary-value problem: (|131|) . y'(0) = 0, and y G L 2 [0,oo). Be- 
cause x I/+2 — > oo as x — > oo, this problem defines a discrete (pure point) spectrum 
S := {Xk\k G Z + } which is real and unbounded, [221 §10.3]. Furthermore, the eigen- 
function associated with A& is, up to a multiplicative constant, 2/2 (• ; A&), and the func- 
tion m has the following spectral resolution: 



-(A) = EA. (135) 



k=l 

-1 



Afc — A 



where cx^ = [J °° ||/2(x; Afc)| 2 <ix] G R. In particular, m is a holomorphic function in 
C — S and A& are the poles of m which are all simple. 19 

Next, consider the following two possibilities: 

1. A + G R or A_ G R: First suppose A + G R, then A_ ^ R and we have 

y_(x; A_) = C(A_)[3/i(x; A_) + m(A_)y 2 (x; A_)], (136) 
where m is given by (|135|) . Eqs. (|129|) . (|13U|) . and (|136|) imply 

y + (0) = y_(0) = C7(A_)m(A_), y' + (0) = -e~ 4 *V(0) = e~ 4 ^(A_), (137) 
and consequently 

y+ (0)-e m »m(\_)y' + (0) = 0. (138) 
In view of (|127|) . which implies A + = e~ 4ie,/ A_, and (|135|) we can express (|138|) as 

Z/+ (0) + x(A + )y , f (0)=0, (139) 

where 

00 

fc=l K 

Next, consider a fixed A + G R. Then because we have the limit point case there 
is at most one linearly independent L 2 -solution y + of 

-^(x)+x^V(x) = A + y + (x). (141) 
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Note that > for all k £ Z + and that S has no accumulation (cluster) point. 



24 



This implies that y*,, which also belongs to L 2 [0, oo) and solves (|141|) . satis- 
fies y+(x)* = e ll y + (x) for some 7 G [0,27r). Inserting this equation in the 
one obtained by taking the complex-conjugate of both sides of (j!39|) and us- 
ing y + (0) ^ ^ y+(0), we have x(KY = In view of d3>> the latter 
relation reads $i(A+) = where 

00 / 1 1 \ 

Hence A + is a real zero of $1. Clearly $1 is a holomorphic function in C — S± \JS^ 
where Sf := {e ±m "\ k \ k G Z+}. Therefore, its zeros (if exist) form a discrete 
set. This in turn means that A + and consequently the eigenvalues E n = e 2l6v \ + 
(associated with this case, if there are any) belong to discrete sets. The same 
argument applies for the case A_ G K. In summary, the eigenvalues that lie on 
the rays: arg(z) = ±2i9 u form a possibly empty discrete subset of C. Next, we 
show that the same holds for the eigenvalues lying outside these rays. 

2. A + ^ E and A_ ^ K: In this case we can use (J132|) to express y± as 

y±(x) = C(A±)[ 2/l (x; A±) + m(A±)y 2 (x; A±)]. (143) 
Substituting this relation in (|129|) and (|13U|) . we obtain 

C(A+)m(A + ) = C(A_)m(A_), C(A_) = -e 4 ^C(A + ). 
These together with (fT27l) . ip35|» . and C(A±) ^ yield 

$ 2 (£„) = e 2ie "m(e 2 ^£ n ) + e^m^" 2 ^ £„) = 0, (144) 

where 

00 / 1 1 \ 

$2(A) ^ ~ £ U-e 2 ^A fc + A-e- 2 ^A fc ) 



a k . (145) 



Therefore the eigenvalues E n are the zeros of $2- 20 Clearly, $ 2 is a holomorphic 
function in C - (S 2 U S 2 + ) where 5^ := {e ±2ie »\ k \k G Z+}. This implies that the 
zeros E n of $2 form a discrete set. Hence the eigenvalues that do not lie on the 
rays arg(z) = ±2i8 u also form a discrete set. 

This completes the proof that the set of all the eigenvalues E n is discrete. □ 



20 Note that in light of l|134|) we have <&2(A*) = $2(A)* = <& 2 (A). Hence the complex-conjugate of every 
zero of $2 is also a zero of $2- This is consistent with the fact that the eigenvalues of H are either real or 
come in complex-conjugate pairs 
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